Energy calibration of the Nal(Tl) calorimeter 
of the SND detector using e+e~ e^e~ events 



M.N.Achasov, ^ D.A.Bukin, T.V.Dimova, V.P.Druzhinin, 
V.B.Golubev, V.N.Ivanchenko, A.A.Korol 

Institute of Nuclear Physics, Novosibirsk, 630090, Russia 



Abstract 

Calibration of the three layer Nal(Tl) spherical calorimeter of the SND detec- 
tor using electron - positron scattering events is described. Energy resolution of 
5%(FWHM/2.36) for 500 MeV photons was achieved. 



1 Introduction 



The SND is a general purpose nonmagnetic detector ( Fig. 1 ) operating at 
VEPP-2M e~^e~ collider in BINP ( Novosibirsk ) in the center of mass energy 
range of 0.2 1.4 GeV. [1,2]. Experimental studies include decays of p,uj,(f) 
mesons and nonresonance hadron production at low energies. Good energy 
resolution for photons in a wide energy range from 30 to 700 MeV is essen- 
tial for suppression of background in reconstruction of intermediate 7r° and 
7] mesons and detection of photons emitted in radiative transitions between 
different quarkonium states. Fast preliminary calibration of the calorimeter 
is based on cosmic muons [3]. It provides reasonable energy resolution of 
5.5%(FWHM/2.36) for 500 MeV photons, but to increase resolution to its 
highest value it is necessary to use experimental events with precisely known 
energies of final particles. In addition, such a process should have clear event 
topology, be well separated from background, and have large cross section. Po- 
tentially suitable processes are e'^e' e+e", e~^e~ — > 77, e~^e~ tt+tt", e+e~ - 
but given the VEPP-2M luminosity of 3 ■ 10^° cm~^s~^ at 1 GeV, only 
e~^e~ e^e" produces enough statistics for calibration of the calorimeter in 
a reasonably short time. 
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2 The SND calorimeter 



SND detector [1] ( Fig. 1) consists of a cylindrical drift chamber, calorimeter, 
and muon system. The thrcc-laycr Nal(Tl) spherical calorimeter, consisting 
of 1632 individual counters ( Fig. 2) is a main part of the detector. It was 
described in [3], so let us mention only some details necessary for description 
of the calibration procedure. 

Calorimeter solid angle coverage in a spherical coordinate system with Z axis 
directed along the electron beam is 18° < 9 < 162° and 0° < < 360°. 
Calorimeter is logically divided into two parts: "small" angles 18° < 6 < 36° 
and 144° <9 < 162°, and "large" angles: 36° < 6' < 144°. The angular dimen- 
sions of crystals at "large" angles are = A6' = 9° and = 18°, A0 = 9° 
at "small" angles. The calorimeter layers arc cmimcrated starting from the 
interaction point - the first one is the nearest to the beam. 

The calorimeter energy resolutions for 500 MeV electrons, photons, and muons 
after primary calibration using cosmic muons are shown in the Table 1 . While 
peak positions in Monte Carlo simulation and experiment agree at a 1% level, 
experimental resolutions are significantly worse than simulated ones. The pos- 
sible explanations of the differences in simulated and experimental spectra 
could be attributed to instabihty of the detector electronics, systematic er- 
rors in the cosmic calibration procedure, and inadequate treatment of nonuni- 
formity of light collection efficiency over the crystal volume in Monte Carlo 
simulation. 

Relative stability of calibration constants in time [3], shows that electronics 
and photodetectors instabilities do not contribute much into experimental 
resolutions. To eliminate systematic biases of cosmic calibration procedure, 
the OFF-LINE calibration based on e+e^ — > e+e^ events was performed. 



3 The calibration method. 

The calibration constants based on e+e~ scattering events could be obtained 
by minimization of the expression: 

F{Q)^Y.ii:uij-Q-Eor, (1) 

j i 

where j - event number, i - crystal number, Uij is an energy deposition in 
ith crystal and jth event, Eq - beam energy, Q - calibration constant for ith 
crystal. Similar calibration procedure was implemented in CLEO II detector 
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[4], but layered structure of the SND calorimeter complicates the task due to 
large fluctuations of energy depositions in the calorimeter layers and energy 
dependence of longitudinal development of electromagnetic shower. Although 
r.m.s. of the energy deposition spectra in an individual crystals are about 
100%, the statistical accuracy of Cj must be high, due to strong correlations 
between energy depositions in different crystals in an event. The drawbacks of 
such a direct method are that the calibration constants are dependent on the 
energy of the electrons and produce biased values of average energy depositions 
in calorimeter layers. 

To avoid such complications the SND calorimeter calibration is based on com- 
parison of coefficients Cf"^ and C^^"^, obtained using simulated and experimen- 
tal electrons of the same energy. 

After cosmic calibration the measured energy deposition in ith crystal t/j can 
be written as Ui = ei-Ei, where Ei is an actual energy deposition and difference 
of ei from unity, characterizes the systematic shift of cosmic calibration for 
ith crystal. To compensate for this shift it is necessary to find corresponding 
correction coefficient Cf^^ — l/e^. These cahbration constants are obtained in 
the following way. First, the functions: 

= E(E CrE^J - E,f and F(Cr ) = Y.(L Cr Uij - Eo)^(2) 

3 i j i 

are minimized over C™'^ and Cf^"^ . Here i is a crystal number, j - event number 
Eij - energy depositions in crystals in Monte Carlo simulation, Uij - measured 
energy depositions in experimental e'^e" events, Eq - beam energy ( the same 
in experimental and simulated events ). The minimums of the functions are 
determined by the following conditions: 

dF{Ci)/dCi = (3) 

Coefficients Cf^^ and C™'^ are the solutions of the set of simultaneous linear 

equations of the form A • c — b, where A is an n x n matrix and its elements 

are — J^EijEmj or J^UijUmj, c and b are vectors of a dimension n with 
j j 

elements d = Cf"" or Cf^ and h = EoJ2 Eij or Eq E Uij. Here n = 1680 is a 

3 3 

total number of crystals in calorimeter, indices i and m - crystals numbers, j - 
the event number. The solutions of the two linear systems satisfy the following 
condition = Ve^, hence Cf = CP/Cf^^ 

The calibration constants Cf'^^ obtained this way have high statistical accuracy, 
are independent of the energy of electrons, and do not produce biases in energy 
depositions in calorimeter layers. 
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For calorimeter calibration e"'"e~ e+e~ events were simulated. In order to 
save simulation time, angular distribution was set to uniform over the solid 
angle. The passage of electrons through the detector was simulated by means 
of the UNIM0D2 code [5]. 

Both experimental and simulated e~^e~ e^e~ events are selected according 
to the same criteria: only two particles must be detected, the total energy 
deposition in the calorimeter is greater than 1.2 • Eq and acoUinearity angle 
is less than 10 degrees. All crystals with energy depositions less than 5 MeV 
are discarded. Remaining crystals are put into calculation of the elements of 
matrix A and vector b. When events processing is finished, the linear system is 
solved using SLAP2 [6] package, i.e. C^^^ and C™'^ are obtained for all crystals 
and then C^^^ are calculated. 



4 Events processing results. 

This calibration procedure was used in the OFF-LINE processing of the data 
collected in 1996 -r- 1997 in the center of mass energy range 0.99 ^ 1.04 GeV 

[1]- 

To obtain C™'^, 50000 simulated events with 500 MeV electrons were processed, 
corresponding to about 150 electrons per crystal. The mean (C™'^) and their 
r.m.s. values cr™^ are listed in the Table 2. In principle, the coefficients may 
depend on the electron energy, layer number, and crystal size. But at this level 
of statistics no significant difference in Cf"^ values for different crystal sizes 
in the "large" angle part is seen. The statistical accuracy cr^*^ of Cf"^ can be 
estimated as cxg'^ = ct'^7(C"^^) ( Table 2 ). 

To obtain C^^^ constants, e'^e~ e+e~ events corresponding to integrated 
luminosity about 130 nb~^ are needed. Such sample contains about 240000 
electrons in the "large" angles part, corresponding to at least 150 electrons per 
crystal. On average, the SND acquires such an integrated luminosity in three 
days of VEPP-2M operation. The mean (C^^), their r.m.s. values and sta- 
tistical accuracy of C^^"^ (c'c ^) are listed in the Table 2 together with the mean 
(C"^^'), their r.m.s. values cr^^^ and statistical errors a(f —\J{cr^c^y + {gq^Y- 

The statistical accuracy of calibration constants for the first two layers is sat- 
isfactory, but for the third layer it is larger than that for cosmic calibration 
[3] . This is due to relatively small energy deposition of electromagnetic show- 
ers in the third layer, combined with large number of hit crystals. Thus the 
infiuence of the accuracy of calibration constants in the third layer on an over- 
all calorimeter resolution for electrons and photons is small. The situation is 
different for muons and charged pions, where relative energy deposition in crys- 
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tals of the third layer is large. In this case high statistical error in calibration 
coefficients significantly increases the widths of energy deposition spectra. So, 
instead of C^^^, the coefficients obtained during cosmic calibration were used 
for the third layer. 

The r.m.s. difference a^^^ in calibration coefficients obtained using cosmic and 

e~^e~ calibration procedures can be estimated as a^^^ = \J [a^^^Y ~ {^CJ^Y- 
This value is about 4% for the first two layers and is less than 5% for the third 
layer. 



5 Energy resolution of the calorimeter. Implementation of the cal- 
ibration procedure 

As a result of e^e~ calibration the calorimeter energy resolutions for 500 MeV 
electrons and photons ( Table 1, Fig. 3 ) were improved by 10%, but still remain 
worse than those expected from Monte Carlo simulation. 

A Monte Carlo simulation of energy deposition was first carried out with a 
uniform description of the nonuniformity of the light collection efficiency over 
the crystal volume. Then, calculations taking into account nonuniform "direct" 
light from a scintillation and diffuse reflection from the crystal boundaries and 
wrapping were performed. The energy deposition U in crystal measured in 
experiment is 

C/(MeV) = C(MeV/pC) -e-v C(l/MeV) • V • E(MeV), (4) 

where E is an energy deposition in the scintillation counter, C - ratio between 
the collected electric charge from photodetector, measured in pC, and energy 
deposition in the units of MeV, v - quantum efficiency of the photodetector 
multiplied by its gain, C, - light yield of the scintillator, ip - light collection 
efficiency, e - electron charge. In general, ip is a. function of coordinates within 
the crystal and depends on reflection coefficient on the crystal boundaries and 
photocathode diameter. 

To take into account the nonuniformity of light collection efficiency, the crys- 
tal response was simulated taking into account a uniform light collection due 
to diffuse reflection and nonuniform "direct" hght collection, depending on 
a solid angle of a photocathode, visible from a scintillation point. The re- 
sults of such simulation the energy distribution width for 500 MeV photons 
is 4.2%(FWIIM/2.36) ( Fig. 4 ). The experimental spectrum is shown in the 
same figure. Calorimeter resolutions for other types of particles are also shown 
in Table 1. The agreement between experiment and simulation became much 
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better and the residual disagreement could be attributed to difference in dif- 
fuse reflection coefficients values in different crystals. 

To study the calorimeter response for photons as a function of photon energy, 
the events e'^e~ — > 77 and e+e~ e~^e~j were used. The kinematic fit of the 
e~^e~ e"'"e~7 events was performed taking into account energy-momentum 
conservation and the reconstructed photon energies E^^^ were obtained. These 
values were compared with direct calorimeter measurements. One could expect 
strong correlation between these values, but study of the simulated e'^e~ — > 
e'^e^j events, where the energy of the photon is precisely known, showed, that 
even for photons in the energy range from 30 up to 150 McV these correlation 
do not change significantly the calorimeter response. The energy range above 
170 MeV was also studied using 27 annihilation events and results agree well 
with those for e'^e~ — > e'^e~7 reaction. The dependence of calorimeter energy 
resolution on photon energies ( Fig. 5 ) was fitted as 

4- 9^ 

aE/E{%) = ■ (5) 

After calorimeter calibration with e^e~ e^e~ events, the photon energies 
turned out to be biased by about 1% ( Table. 1, Fig. 4 ). In oder to compensate 
this bias, the calibration coefficients for photons were corrected accordingly. 
The distributions over two-photon invariant masses m-y-y in ^ 777 and Ks — > 
7r°7r° decays after such correction are shown at Fig. 6 and 7. Peaks at 7r° and 
r] mesons masses are clearly seen. 

Relative drift of cahbration coefficients is shown in Fig. 8. It can be seen, 
that for a time period between consecutive calibrations, the mean shift of the 
coefficients is about 1% and r.m.s. of their random spread is about 2.5% for 
the first two layers and about 5% for the third layer. 



6 Conclusion 

Using the described procedure for the SND calorimeter calibration the sta- 
tistical accuracy of 2% in calibration constants for the first two layers was 
achieved. The final resolution for photons varies from 10% at 50 MeV to 5% 
at 500 MeV. 
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Table 1 

Calorimeter response for 500 MeV electrons, photons and muons. Eq = 500 MeV — 
energy of the particles, E'^'f''* — measured energies in the calorimeter for electrons, 
photons and energy deposition for muons respectively, a — FWHM/2.36 of distri- 
bution over E/Eq, Peak — the peak position in the distribution over E/Eq. EXPl 
and EXP2 — experimental distributions after cosmic and e~^e~ calibrations re- 
spectively. MCI and MC2 — distributions in Monte Carlo simulation without and 
with nonuniformity of light collection over crystal volume was taken into account 
respectively. 





EXPl 


EXP2 


MCI 


MC2 


Peak a{%) 


Peak a{%) 


Peak a{%) 


Peak a{%) 


E^/Eo 


0.99 5.4 


1 4.7 


0.99 3.5 


1 4.2 


EVEo 


1 5.4 


1.01 5.0 


1 3.7 


0.99 4.2 


Ei'/Eo 


0.33 8 


0.33 8 


0.33 5 


0.34 7 



Table 2 



(jmc^ (jcxp^ (jca\ coefficients and their statistical accuracy in calorimeter layers for 
crystals at the "large" angle zone 



layer number 














I 


1.04 0.02 


1.9 


1.07 0.04 


1.5 


1.03 0.04 


2.3 


II 


1.05 0.02 


1.8 


1.05 0.04 


1.5 


1 0.04 


2.3 


III 


1.29 0.06 


5 


1.27 0.07 


4 


0.98 0.08 


6.4 
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20 40 60 80 100 cm 

Fig. 1. SND detector, section along the beams; 1 — beam pipe, 2 — drift chambers 
3 — scintillation counters, 4 — light guides, 5 — PMTs, 6 — Nal(Tl) crystals, 7 — 
vacuum phototriodes, 8 — iron absorber, 9 — streamer tubes, 10 — 1cm iron plates 
11 — scintillation counters, 12 and 13 — elements of collider magnetic system. 
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Fig. 3. Energy spectra for 500 MeV electrons; Eq = 500 MeV - beam energy, E - 
measured energy. EXP2 - distribution after e^e^ calibration, EXPl - distribution 
after cosmic calibration. 
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Fig. 5. Dependence of the calorimeter energy resolution on the photon energy. 
E - photon energy, ge/E - energy resolution of the calorimeter obtained using 
e^e~ — > 77 ( dots ) and e+e~ e^e~^ ( circles )reactions. The error bars show 
only statistical errors. 




Fig. 6. Two photon invariant mass distribution in the experimental — > KsKl, 
Kg 7r°7r° events. Line - asymmetric Gaussian fit. 




Fig. 7. Two photon invariant mass distribution in experimental (j) ^ rj^ events. Line 
- asymmetric Gaussian fit. 
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Fig. 8. The calibration coefficients spread. Points — average ratio of a current cali- 
bration result to the preceding one, error bars — FWHM/2.36 of the distributions 

of these ratios over the whole calorimeter layer. Horizontal axis shows the time 
elapsed from the first calibration. Shown are the results for the second calorimeter 
layer. Other layers behave similarly. 
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